Enhancement of persistent currents due to confinement in metallic samples 
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Confinement and surface roughness (SR) effects on the magnitude of the persistent current are anal- 
ized for ballistic bidimensional metallic samples. Depending on the particular geometry, localized 
border states can show up at half filling. These border states contribute coherently to the persistent 
current and its magnitude is enhanced with respect to their value in the absent of confinement. A 
linear scaling of the typical current It yp with the number of conduction channels M is obtained. 
This result is robust with respect to changes in the relevant lengths of the samples and to the SR. 
Possible links of our results to experiments are also discussed. 



PACS 72.23.+R 

In recent years, advances in nanotechnology made pos- 
sible to design mesoscopic samples in which the carriers 
are confined and mainly scattered by the boundaries of 
the system Q. In this situation, where the elastic mean 
free path is much larger than the system size, the sam- 
ple can be considered in the ballistic regime in order to 
compute the relevant transport quantities. In addition, 
the potential that confines the carriers into the meso- 
scopic device is not perfectly controlled in experiments. 
Therefore, even when the bulk disorder is absent, surface 
roughness (SR) is present and as a consequence the carri- 
ers are irregularly scattered at the borders of the sample. 

When the sample is threaded by a magnetic flux per- 
sistent currents originate |^[|. Persistent currents (PC) 
have been studied in confined bidimensional cylin- 
drical systems without SR employing discrete models 
with periodic boundary conditions in one of the direc- 
tions. 

Without reference to PC, a discrete lattice has been 
considered recently in order to take into account SR in 
metallic systems J7j. In the latter model, a kind of su- 
perficial disorder was introduced through the absence of 
hopping to and from L sites chosen at random among the 
L 2 sites defining the lattice. 

In this letter we study the PC for a model in which 
the confinement and the SR effects can be simultaneously 
analized. We consider circular (centred and non-centred) 
and square bidimensional loops on a discrete square lat- 
tice. Therefore the confinement breaks completely the 
translational symmetry of the underlying lattice. This 
introduces a qualitative difference with respect to cylin- 
drically shaped systems, in which the translational sym- 
metry is preserved in the azimuthal direction and the two 
degrees of freedom arc decoupled. In addition, the sam- 
ples studied by us have the geometry of those considered 
in the experiment of Ref. || . 

Although much work has been done concerning the in- 
fluence of bulk disorder on persistent currents j^] , less is 
known about the ballistic regime. Our goal is to char- 
acterize the effect of the confinement and the SR on the 



magnitude of the persistent current in clean metallic sam- 
ples. Confinement effects in the context of orbital mag- 
netism were previously considered for semiconductor de- 
vices Jl§. 

The total current can be calculated at zero tempera- 
ture as / ~ dE/dcf), with E the total energy of the system 

and 



the magnetic flux. The value of I typ = y J I 2 d(j> 
depends on the number of open channels present in the 
system, M. For free electron models with unbounded en- 
ergy spectrum it is known that I r 



M |||10|. This 
behavior can be understood taking into account the prop- 
erties of the crossings (or quasicrossings) that appear in 
the spectrum when the flux is varied and frustrate the co- 
herent increase of the current . As expected, the same 
scaling with M holds for tight binding models away from 
half filling || . At half filling and for cylindrically shaped 
systems, it was found that very particular geometrical 
configurations of the sample could give rise to coherent 
contributions of all the channels, but those sharp config- 
urations would be very difficult to obtain in an experi- 
mental set up (HJ6). 

We consider a N x N cluster of the square lattice. Tak- 
ing a single atomic level per lattice site the Hamiltonian 
is, H = Em e m|m >< m| + E m; k*m,k| m >< k| , 
where (to, n) = m labels the coordinates of the sites 
in the lattice. All the lengths are given in units of the 
lattice constant a. The hopping integrals i m;m ' are re- 
stricted to nearest neighbours. Assuming that the vec- 
tor potential has only azimuthal component we take 
t m:m+ i = t exp (i J i ^ 1+1 A.dl), where the phase is mea- 
sured in units of the quantum flux <j> and 1 is a vector 
that points from the site m to any of its four nearest 
neighbours. With the on-site energies, e m , we simulate 
the confining potential in order to get the required profile 
of the sample. We will first focus on the circular centred 
geometry. As the confining potential and the underlying 
lattice have different symmetry, SR is present in these 
samples but the rotational symmetry in 7r/2 is preserved. 
The lattice parameter a represents a cut off for the al- 
lowed wavelengths and therefore the SR will be relevant 
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FIG. 1. (a) Energy levels at the bottom of the spectrum 
as a function of the rescaled magnetic flux a for sample SI 
together with a typical charge distribution, (b) Idem as (a) 
but for the quarter filling region of the spectrum. 



near half filling, that is in the metallic regime. The two 
length scales are chosen as the internal radius r and the 
width of the sample W . While the second determines 
M, the ratio s = W/r controls the strength of SR. For 
s » 1 and for carriers with a wavelength of the order of 
a, two succesive scattering events by the inner and outer 
boundaries of the sample are almost uncorrelated. On 
the other hand, for s « 1 those events are strongly cor- 
related. Therefore the eigenstates corresponding to these 
particular wavelengths will evolve from generic extended 
ones to strongly localized as the value of s is reduced. 

In Fig. [I] we show for sample 51 (with W — 40 and 
r = 20) the bottom (Fig. |l](a)) and the quarter filling 
(Fig. 0(b)) regions of the energy spectrum as a function 
of the normalized flux a = 4>/4>q together with a generic 
charge distribution for each region. Figure |l|(a) corre- 
sponds to the largest wavelengths, that are the less sen- 
sitive to the effect of the SR. This region of the spectrum 
looks qualitatively similar to the spectrum of an inte- 
grable Aharonov-Bohm annular billiard |Tl|| . Moreover, 
the generic charge distribution corresponds to eigenstates 
with well defined value of the angular quantum number. 

In the quarter filling region the typical wavelengths are 
larger than a but small enough to make the states sensi- 
tive to an effective polygonal boundary. This implies the 
lost of the continuous symmetry characteristic of the very 
large wavelengths. Therefore, as the energy rises towards 
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FIG. 2. hyp as a function of N p for sample 51. Inset: 
Idem as (a) for a bidimensional square loop with outer square 
boundary of 20X20 and inner one of 6X6. Since the mean 
perimeter of the square loop doubles the mean perimeter of 
51 the maximum value of I typ in the inset is almost twice the 
value for 51. 



half filling, avoided crossings show up. The correspond- 
ing charge distributions arc quite similar to those for a 
polygonal billiard and are extended states in general (see 
Fig. 1(b)). 

In any case, when an occuppied level crosses (or qua- 
sicrosses) an empty one, the total current / exhibits a 
discontinuity (or abrupt oscillation) ||. As it is shown 
in Fig. ||, below half filling, I typ is a highly fluctuating 
function of the number of particles, N p . This is a conse- 
quence of the random distribution of crossings (or quasi- 
crossings) that appears in these regions of the spectrum 
|nj. The upper region of the spectrum, close to half 
filling, is shown in Fig. |3|(a) for 51. In this region of en- 
ergies the wavelengths are of the order of a. A bunch of 
quasi-degenerate states appears at zero flux at half fill- 
ing which manifestes as a peak in the density of states 
(DOS) at zero energy. For other circular centred samples 
a qualitatively similar behavior was found. In a bidimen- 
sional sample, a peak in the DOS at zero energy could 
appear for a very specific relation between the two rele- 
vant lengths of the sample. For example, for a cylinder 
of circumference L (in units of a) and with M trans- 
verse channels the condition is L = 2 (p/j) (M + 1) [Q, 
where p and j are coprime integers. For a rectangular 
sample, i. e with Dirichlet boundary conditions in both di- 
rections, with L x K sites the corresponding formula is 
L + 1 = 2 (p/j) (K + 1) . As a consequence, for a square 
shaped sample the above relation is always satisfied and 
a peak in the DOS manifestes. We expect this result 
to remain valid in general for the samples considered in 
this work, in which Dirichlet boundary conditions are 
imposed on the inner and outer borders, providing the 
confinement preserves the rotational symmetry in tt/2. 
This is confirmed by the numerical results presented in 
Figs. ^ for the square loops. Therefore, independently of 
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FIG. 3. (a) Energy levels close to half filling as a function 
of the rescaled magnetic flux a for 51. (b) Typical charge dis- 
tributions for angularly localized states (left side) and radially 
localized states (right side). 



the number of channels in the loop, the square and cen- 
tred circular geometries have a peak in the DOS at zero 
energy. For the circular loops, as a consequence of the 
SR, the states that conform the peak in the DOS have 
different characteristics. Some of them are angularly lo- 
calized states (see the left side of Fig. |](b)). They are not 
sensitive to the magnetic flux and correspond to the flat 
lines in the spectrum. The other type of states are mostly 
radially localized and look like whispering gallery modes. 
These states are originated by the correlated scattering 
with the outer and inner boundary of the sample that al- 
lows that the charge distribution does not fill completely 
the loop and the state remains localized. Althoug local- 
ized in one direction, these states go around the sample 
being sensitive to the magnetic flux and carrying a finite 
current. Moreover, they have a quasi well defined value 
of the angular quantum number that characterizes the 
slope of the eigenenergies as a function of the flux. The 
right side of fig. [|(b) displays the characteristic charge 
distributions for these states. 

Opposite to what happens for systems with diagonal 
disorder, in the present model the localization is obtained 
without introducing an additional scale of energy. Only 
remains the scale defined by the kinetic energy. At zero 
flux, the energy of the localized states will be very close 
to zero as a consequence of the symmetry and the finite 
character of the spectrum. Therefore, the existence of lo- 
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FIG. 4. DOS for square loops with different widths W. 
Notice the peak presents in the DOS at half filling. Inset: 
DOS for the rectangular samples. There is no bunching of 
states at zero energy in general. 



calized states in the system is manifested as a peak in the 
DOS that is equivalent to a finite gap in the spectrum of 
eigenergies. Before estimating the magnitude of the gap 
we remark that since the eigenenergies are degenerate at 
a = but with quite different values of their slopes, each 
of the radially localized states must belong to a differ- 
ent conduction channel. Being radially localized, these 
states are quasi one dimensional and the energy gap is 
estimated as A g = (dE/dk) Afc > (7 A)^ 1 , where Afc is 
the inverse of the mean radius of the ring and 7 is a finite 
fraction of W which measures the degree of spatial local- 
ization of the state. In order to check the above result we 
have done numerical calculations for samples with same 
area A, obtaining the same value of A g . The gap has 
an important consequence on the values of It yp close to 
half filling. As it is shown in Fig. @, I typ increases mono- 
tonically with the number of occupied levels (number of 
particles) N p . 

Some comments must be made in connection to the 
non-centred circular loops. Although in this case the SR 
does not preserve the symmetry in 7r/2, it affects the in- 
ner and outer boundary of the samples in the same way. 
This allows the correlated scattering between the inner 
and outer borders and therefore angularly and radially 
localized states also raise up at half filling, as it is il- 
lustrated in the spectrum of Fig. |j| In the case of the 
square loops, although the results are qualitatively simi- 
lar to those obtained for the circular samples, some im- 
portant differences must be noticed: i) The square loops 
do not have SR and therefore there are not avoided cross- 
ings in the spectrum as a function of the flux, ii) Border 
states carrying a finite current were also found (analogous 
to the radially localized states of the circular samples). 
Therefore It yp increases coherently near half filling as it 
is shown in the iinset of Fig. ^|. Hi) As a consequence of 
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FIG. 5. (a) DOS for a non-centred circular loop of W = 9. 
(b) Energy espectrum as function of the flux for the same 
sample. Flat lines at zero energy are associated with angu- 
lary localized states. Flux depending eigenvalues begining at 
zero energy are associated with radially localized states. The 
number of the latter (former) increases (decreases) with W. 

the absence of SR, angularly localized states do not exist 
at half filling (notice the absence of the plateau in I typ 
at half filling) . 

For rectangular samples the localized border states 
could occur for very specific relations between it lengths, 
and in general the spectrum has not bunching of states 
at zero energy. In other words, a peak in the DOS is 
not expected at half filling for the rectangular geometry, 
as it is shown in the inset of Fig. |^. We also find, in 
general, that when the localized border states that carry 
a finite current become dominant over the angularly lo- 
calized states the scaling of It yp with M at half filling is 
substantially enhanced with respect to the known scal- 
ing law out of half filling {Ityp ~ V~M)- Fig. || shows the 
actual scaling law obtained for the square loops. Notice 
that the linear increase of It yp with M is consistent with 
the fact that each of the localized states bunched at zero 
energy belongs to a different conduction channel. 

The effect described in the present work can help to un- 
derstand some puzzles about PC in metallic rings. This 
is the case of the experiment of Ref. [|| in which the 
PC of isolated Au rings was measured to be a factor 30 
to 100 times larger than the theoretical estimates. In 
addition, in that experiment the reported value of the 
signal for the circular loop was 5 times larger than for 
the rectangular loop (although the predicted theoretical 
values were roughly the same). Our results, even with- 
out providing a definitive explanation of the discrepan- 
cies between theory and experiment, are in line with the 
strong enhancement of the current reported in it and pro- 
vides the way to explain the observed differences between 
the circular and the rectangular samples. In fact, as we 
previously described, for the circular loop we predict a 
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FIG. 6. Linear scaling of It yp with the number of chan- 
nels M, at half filling, for square loops with the same mean 
perimeter. 

peak in the DOS and therefore the enhancement of the 
PC. On the other hand, as we already discussed, for the 
rectangular sample the enhancement is not expected in 
general. As a last point, we would like to remark that 
the enhancement of the current predicted in this work 
does not have a critical dependence on the dimensions 
and on the SR of the samples, making their experimen- 
tal realization feasible. Unfortunately we are not aware 
of new measurements with other sample sizes in order to 
confirm or not our predictions. 
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